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SUMMARY 

A solution technique has been developed for solving the multiple-blade- 
element, surface-of-revolutlon, blade-to-blade flow problem In turbomachinery. 
The calculation solves approximate flow equations which Include the effects of 
compressibility, radius change, blade-row rotation, and variable stream sheet 
thickness. An Integral equation solution (l.e., panel method) Is used to 
solve the equations. A description of the computer code and computer code 
Input Is given In this report. 


INTRODUCTION 

A computer code has been written to analyze flow on turbomachinery blade- 
to-blade surfaces. The code has the capability of analyzing blade rows which 
are made up of several different blade shapes and/or spaclngs. This capability 
makes the code useful for calculating flows In centrifugal machinery with 
splitter blades or In mls-tuned blade rows, where the blade spacing varies. 

The number of varying blade shapes that can be used In a blade-row design Is 
not limited by the code. However, the amount of computer storage available 
does Impose a practical limit on any problem. The code being described here Is 
dimensioned to handle from one to four bodies In cascade. After redimensioning, 
this same code has been used to calculate a flow problem with as many as 15 
bodies. 

The basic method used to solve the flow field was described In reference 1. 
However, the method of reference 1 was extended In the code reported here to 
Include multiple-body cascades. Briefly, the solution method solves approxi- 
mate governing equations for the blade-to-blade, steady-state flow problem In 
turbomachinery. The effects of compressibility, radius change, blade-row rota- 
tion, and variable stream sheet thickness are Included. The working fluid Is 
assumed to be Invlscld, Irrotatlonal , and a perfect gas. Comparisons with other 
solution techniques and experimental data Is given In reference 1. Overall, 
the solution produces good results for subsonic flow, but the accuracy decrea- 
ses for high subsonic and supersonic flows. 


CODE DESCRIPTION 
General 

The multiple-body-panel code Is written primarily In FORTRAN IV, and uses 
double precision arithmetic. Storage for the four-body version of the code. 
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which accompanies this report, Is approximately 650 000 words when run In 
double precision. Typical run times for the code on the IBM 370/3033 computer 
at NASA Lewis Is 25 CPU sec. A slightly different version of the code which 
runs on the CRAY IS computer at NASA Lewis requires less than 5 CPU sec/run. 

The computer code accompanying this report Is essentially the version of 
the code that runs on the CRAY computer. However, coding that Is CRAY computer 
dependent and allows vectorlzatlon of the calculation has been replaced with 
scalar coding which Is machine Independent. Also, calls to the Society of 
Industrial and Applied Mathematics (SIAM) matrix solvers contained In their 
LINPACK software package have been replaced by calls to an In-house generated 
solver which Is Included with the code. The result of these code changes Is 
a code which Is slower running than the version of the code used at NASA 
Lewis, but Is more machine and system Independent. If the user has access to 
SIAM'S LINPACK software, Instructions on how to Impllment It Into the code are 
given In the comment statements of subroutine PLMS0L. 

Some machine dependency was left In the code. The code Is set up with the 
coding for variable storage arrays and dimensions separated from the main body 
of the code. On both the IBM and CRAY computer systems at NASA Lewis, this 
coding Is Inserted Into the main code by the compiler. This compiling option 
allows the user to adjust the storage requirements of the code and to Increase 
(or decrease) the number of bodies which can be considered. The option Is 
available on many computing systems, and so the code's structure was left In 
this form. If the user's system does not have this capability, he can physi- 
cally dimension the arrays and Insert the coding Into the program before com- 
piling. 


Code Structure 

The code on the tape as given to COSMIC consists of several short pieces 
of coding followed by the main code and Its accompanying subroutines. The 
short pieces of coding contain the variable array and dimensioning Information. 
These segments are Inserted Into the main coding during compilation. This 
coding appears first on the tape. Each segment Is preceded by a line that 
reads "*CDK" followed by the name of the segment. Five of these segments are 
used In the code. Each segment Is Inserted Into the program wherever the cod- 
ing "*CALL“ followed by the name of the segment appears. The main program and 
subroutines are given next on the tape. A line reading "*DK" followed by the 
name of the routine precedes each subroutine and the main program. The multi- 
ple body panel code consists of the main program plus 24 subroutines. 

The multiple body panel code Is written In a modular style. No calcula- 
tions or Input/output takes place In the main program. The main program only 
serves as a calling routine. Each subroutine then performs a particular func- 
tion In the overall solution. The various parts of the code are listed below 
along with a brief description of the coding's function. Items are listed In 
the order that they should appear on the tape. 
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Variable Arrays and Dimension Information Coding (Code to be Inserted Into 

the Main Program Before Compiling) 

BASBLK - Contains all geometry and flow solution variable common blocks. 

Also contains parameter statements which establishes the size of 
the arrays. 

DIMEN - Sets values of specific array dimensions which are used In 
read/write statements and are passed to variable dimension 
subroutines. 

REG2 - Contains arrays to be used In the calculation of velocities In 
the flow field (off-body velocity calculations). 

REG3 - Contains bulk storage arrays used to restore some of the 

two-dimensional arrays of BASBLK as one-dimensional arrays for 
calculations. Also contains arrays used for temporary storage of 
Intermediate results. 

REG5 - Contains temporary storage arrays for use In subroutine PLMSOL. 


Main Program and Subroutine Coding 

MLTPNL - Main program acts as a calling routine for subroutines. 

PLMIPT - Reads basic Input data file and writes labeled echo check. 

Called by MLTPNL. 

BLADE - Generates blade shapes by patching together spline curves (Input 
option IBLD =1). Called by PLMIPT. 

PLMTRF - Performs the geometric transformation from a surface of 
revolution to a plane cascade. Called by MLTPNL. 

ALPOCL - Calculates the tralllng-edge-statlon flow angle from the 

downstream flow determined during the Kutta condition Iteration. 
Called by PLMSOL. 

PLMFRT - Partially transforms the solution velocity back to physical 
velocity. Called by PLMSOL. 

GEOOPT - Outputs the transformed and physical geometry of the bodies to 

both the output listing file and the plot file. Called by MLTPNL 

CMPOPT - Completes transformation of the solution velocity back to 

physical velocity. Makes additional compresslbllty corrections 
(Input option ISL =1). Outputs final blade-surface flow 
solution to both output listing and plot files. Called by PLMSOL 
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RBINT - Calculates the one-dimensional, compressible flow used In approxi- 
mate governing equations. Also performs Integrations to get con- 
tributions of the compressibility, blade-row rotation, and stream 
sheet thickness variation to the flow solution. Called by PLMSOL. 

PLMSOL - Sets up the system of equations to be solved. Acts as a secondary 
calling routine which brings together all the subroutines and 
Information used to solve the matrix of equations for the flow. 
Called by MLTPNL. 

PLMGEO - Calculates geometric parameters and variables used In the 
solution. Called by MLTPNL. 

KCGEOM - Calculates the points downstream of the trailing edges and the 

flow angles to be used In the Kutta condition (Input option ICSCD 
= 0, 2, or 3). Called by PLMGEO. 

VORTEX - Calculates Influence of a linear vortex singularity distribution. 
Called by ABCOEF and OFFBOD. 

ABCOEF - Acts as a secondary calling routine to call the subroutines which 
calculate the Influence coefficients of the Integral-equation 
solution. Called by MLTPNL. 

JUMP - Calculates the Influence coefficients for a discontinuous change 
In a vortex singularity distribution which Is due to discontin- 
uous changes In the slope of the surface of a body. Called by 
ABCOEF and OFFBOD. 

SOURCE - Calculates Influence coefficients of a linear-source singularity 
distribution. Called by ABCOEF and OFFBOD. 

VELCAL - Calculates the velocity Influence coefficient due to source and 
vortex distributions for Isolated bodies. Called by ABCOEF and 
OFFBOD. 

CASCAD - Calculates the velocity Influence coefficients due to source and 
vortex distributions for a cascade of bodies. Called by ABCOEF 
and OFFBOD. 

BNDCND - Calculates the far-stream boundary conditions from the Input flow 
conditions at the blade leading- and tralllng-edge stations. 

Called by PLMSOL. 

POTOPT - Outputs the transformed velocities and singularity strength 

distribution. Used primarily for solution debugging. Called by 
PLMSOL. 

SPLINT - Calculates a spline curve fit and Interpolation from a given 

array of data points. Assumes the second derivative of the curve 
goes to zero at the end points. Also has an entry point SPLENT 
for multiple Interpolations along the same curve. Called by 
BLADE, PLMGEO, PLMTRF, PLMSOL, BNDCND, RBINT, and OFFBOD. 
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SPINSL - Calculates a spline curve fit and Interpolation from a given array 
of data points. Uses the given slope at the curve ends to calcu- 
late the spline curve fit. Called by BLADE. 

OFFBOD - Reads coordinates of points that the user has specified. Calcu- 
lates the flow velocity and direction at those points. Writes 
out body velocity Information to both the output listing and plot 
files. Called by CMPOPT. 

CHLSKY - Solves a set of linear equations by upper and lower triangular 

matrix decomposition. Has an entry point called CH6CNT for repeat 
solutions of the same system of equations with a different right 
hand side. Called by PLMSOL. 

CFCALC - Calculates the tangential and normal force coefficients for each 
body. Also calculates lift, drag, and moment coefficients. 

Outputs results to the output listing file. Called by CMPOPT. 


Input, Output, and Storage Files 

The multiple-body-panel code Is reasonably straightforward to set up and 
run. The Input data file Is column formated. All numerical data Is either 
8F10.0 for floating point or 1615 for fixed point data. Floating and fixed 
point data are not mixed on any Input line. Details of the Input data file 
are given In a separate section of this report. Input Is read In from unit 
05, and printed output Is sent to unit 06. In addition, an unformated binary 
output data set can be sent to unit 60. This binary data set contains geo- 
metric and flow solution Information. It Is Intended for use In any graphics 
code that the user may develop. A description of the Information In this 
"graphics" file will be given later. In order to reduce the In-core storage 
needed to run the code, extensive use Is made of mass storage devices. This 
requires that units 10 to 20 be assigned when running the code. These units 
are used strictly as scratch storage and should be released upon completion of 
the calculation. 

The use of mass storage (units 10 to 20) Is a trade-off of run time for 
computer storage. If the user wants a faster running code and has a large- 
memory computer system available, he may want to eliminate the out-of-core 
reads and writes. The experience of the author Is that this Increases the 
required In-core storage by about a factor of five for the four-body problem 
(from 0.65 to 4 M words). 

The conversion of the code to run using only In-core storage Is not a 
simple task. First, the user must locate all the Initial writes to units 10 
to 20. These are located In subroutines ABC0EF, OFFBOD, and PLMSOL. Second, 
the user must define two-dimensional arrays to store the values that had prev- 
iously been written out and read In from out-of-core storage. Third, the user 
must substitute the new variables Into the solution equations with the proper 
Indices. Lastly, the user must remove all read and write statements to units 
10 to 20. 
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Code-Generated Diagnostic Messages 

Several messages can appear In the printed output from the multiple-body 
code. The name of the subroutine generating the message Is appended to the 
end of the message. The messages are listed together along with a discussion 
of what occurs In the calculation when the message appears. 

1. THE AVERAGE PANEL IS SMALLER THAN THE ADJACENT L. E. OR T. E. PANELS - NS 
OR NE IS LESS THAN ZERO. THE PERCENTAGE OF POINTS IN THE LEADING EDGE HAS BEEN 
INCREASED. PLE IS NOW - BLADE 

The blade shape generator was set up assuming that the elements used to 
describe the leading and trailing edges would be smaller than the elements 
used to describe the blade surface. If the user specifies Input parameters 
that violate this condition, the code will Increase the number of elements 
used to describe the leading edge until the leading-edge elements are smaller 
than the surface elements. Everytlme the message Is printed, 5 percent of the 
total number of elements to be used In describing the blade are added to the 
leading-edge elements. Increasing the number of leading-edge points does not 
always allow the solution to continue. If the tralllng-edge radius Is large 
compared to the leading edge and has panels that are larger than the average 
blade surface panel, then the blade shape generator will fall to give a distri- 
bution of panel elements that will allow the calculation to proceed. In these 
cases, the user will have to provide the code with a descrete point description 
of the Input body. 

2. THE DOWNSTREAM FLOW CALCULATION IN SUBROUTINE ALPOCL IS IN ERROR - PLMSOL 

This error message can appear In solutions using the Kutta condition option. 
It Indicates that the downstream flow angle, calculated using the blade-row 
circulation, exceeds the maximum possible turning angle for the flow. The 
downstream angle Is set to the maximum value so that the calculation may 
proceed. 

3. KUTTA CONDITION HAS FAILED TO FIND A SOLUTION IN 100 ITERATIONS. ERROR IN 
THE INLET FLOW ANGLE IS _ RADIANS - PLMSOL 

The Kutta condition Iteration (ICSCD = 0, 2, or 3 In the Input) has failed 
to converge to the given Input flow. The calculation will continue using the 
flow field at the 100th Iteration. If the error Is large, the calculated 
results will be Invalid. 

4. ITERATION FOR (UPSTREAM/OUTLET/DOWNSTREAM) FLOW CONDITIONS HAS FAILED TO 
CONVERGE. ERROR IS - BNDCND 

Using the given Inlet flow conditions, the solution calculates a mass flow 
for the blade passage. Flow conditions at the upstream, downstream, and outlet 
blade-row stations are then calculated. If the Iteration procedure Involved 
In solving the equations does not converge In 100 Iterations, the above message 
Is printed with the appropriate station noted, and the calculation Is allowed 
to continue. 
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5. SPECIFIED INLET FLOW EXCEEDS CHOKING MASS FLOW - BNDCND 

If the given Inlet flow conditions define a flow that exceeds the choking 
mass flow, the above message Is printed. The calculation continues but the 
solution Is Invalid. 

6. SPECIFIED OUTLET FLOW EXCEEDS CHOKING FLOW - BNDCND 

7. GIVEN OUTLET FLOW ANGLE EXCEEDS MAXIMUM ALLOWABLE - IT HAS BEEN LOWERED TO 
- BNDCND 

The above two messages are printed when the outlet flow angle produces a 
flow that Is greater then the maximum flow possible. The flow angle Is adjus- 
ted to produce the maximum flow at the exit, and the calculation continues. 

The messages can be disregarded for Kutta condition calculations, since the 
final outlet flow conditions are determined later In the calculation. 

8. MAXIMUM OR NEGATIVE VELOCITY ENCOUNTERED AT (UPSTREAM/DOWNSTREAM) BOUNDARY 
- (UPSTREAM/DOWNSTREAM) RADIUS TOO SMALL ? - BNDCND 

This message Indicates that the critical velocity ratio has exceeded the 
maximum value at either of the far-stream stations as noted. This can occur 
In radial machinery If the far-stream boundary condition Is taken at too small 
a radius. For rotors, the error can also occur If the radius Is too large. 

The solution continues, but It Is Invalid. The user can avoid the error by 
moving the boundary station (MRSP(l) or MRSP(NMN) In the Input) In closer to 
the blade row. 

9. AVERAGE MASS FLOW CALCULATION EXCEEDS CHOKING FLOW AT M = . CALCULATED 
AVERAGE MASS FLOW IS TIMES THE CHOKING MASS FLOW. THE FLOW TURNING ANGLE 
HAS BEEN REDUCED TO ALLOW THE CHOKING MASS FLOW TO PASS - RBINT 

This message can occur for calculations using the blade-passage geometry 
and one-dimensional continuity to calculate the mean flow estimates. The cal- 
culation Is allowed to continue, but the mean-flow turning angle Is set to the 
maximum turning value Instead of using the average blade angle. 

10. AVERAGE MASS FLOW HAS FAILED TO CONVERGE. THE REMAINING ERROR IS . - 
RBINT 

This message can occur for calculations of the mean flow estimate from the 
one-dimensional continuity and the blade-passage geometry. It Indicates that 
the Iteration procedure used to calculate the velocity and density has failed 
to converge In 100 Iterations. The calculation Is allowed to continue without 
further Iteration. 

11. (UPSTREAM/INLET/OUTLET/DOWNSTREAM) VALUE OF THE W/WCR ARRAY IS INCONSISTENT 
WITH VELOCITY DIAGRAM - RBINT 

This message can appear In calculations when the user has given the crit- 
ical velocity ratio distribution through the blade row. It warns the user 
that at the location noted In the error message the given mean-flow critical 
velocity distribution differs from the critical-velocity-ratio values calcu- 
lated by the code at the boundary-condition stations by more than 0.01. The 
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calculation continues. The given critical-velocity-ratio distribution Is used 
to estimate the flow compressibility, and the calculated values are used as 
the boundary conditions for the solution. 

12. MASS FLOW CALCULATION AT THE THROAT EXCEEDS THE CHOKING VALUE - RBINT 

This message only appears for a single body In cascade problem when the 
user has set IMN = 2 or 3. The mean flow Is set equal to the choking value at 
the code-calculated throat and the calculation continues. 


INPUT FILE DESCRIPTION 

The Input data file for the multiple-body-panel code consists of two parts. 
The basic Input data describes the problem geometry and flow conditions. If 
only the basic Input data Is given, the flow velocity Is calculated only on 
the body surfaces. The second part of the Input file Is blocks of data which 
contain coordinates of points In the flow field where the user wants the flow 
velocity and direction to be calculated. This off-body velocity-calculation 
Input-data file Is added to the basic Input-data file when velocities In the 
flow field are needed. 


Basic Input File 

The basic Input-data file Is described In this section. The actual form 
of the Input data varies according to the Input options selected by the user. 

Item 1 - Data set title - up to 80 characters. (Format 20A4) 

Item 2 - Calculation control parameters. (Format 1615) 

NBODY NCS ICSCD IMN NMN ISL 

NBODY - Number of separate bodies or blade shapes which make up a cascade 
group. Minimum number 1, maximum number 4. Code can be redimen- 
sioned to allow more or fewer bodies. 


NCS - Number of flow variation cases to be run. Maximum number 8. 


ICSCD - Type of calculation to be made. 


ICSCD = 0, 


Isolated multi-element airfoil. Kutta condition used to set 
circulation for each element. 


ICSCD = 1, Multi-element cascade solution outlet flow angle given. Circu- 
lation split for each body element Is specified In Item 7. 

ICSCD = 2, Multi-element cascade solution Kutta condition used to split 
the circulation between the cascade body elements. Tralllng- 
edge flow angle taken to be the tralllng-edge angle bisector. 

ICSCD = 3, Same as ICSCD = 2, except the user specifies In Item 7 the 
tralllng-edge flow angle to be used in the Kutta condition. 


8 



IMN - Type of flow path description and mean flow calculation. 

IMN = 0, Radius and stream sheet thickness are given as functions of the 
meridional coordinate. Mean flow Is calculated using given blade 
geometry. Only subsonic flows (with local supersonic pockets) 
can be calculated. 

IMN = 1 , In addition to radius and stream sheet thickness, the mean 

critical velocity ratio Is also given as a function of meridional 
coordinates. All flow ranges. Including supersonic flows, can be 
calculated. 

IMN = 2, For single body In cascade only (NBODY = 1, ICSCD = 1). The 

radii and stream sheet thicknesses are given. Program finds the 
blade passage throat and calculates mean velocity at the throat. 
If the Inlet critical velocity ratio Is greater than the outlet 
critical velocity ratio, a linear distribution Is calculated from 
the Inlet to* the throat, and the blade geometry Is used to calcu- 
late the mean flow from the throat to the outlet. If the Inlet 
critical velocity ratio Is less than the outlet critical velocity 
ratio, the blade geometry Is used to calculate the mean flow from 
the Inlet to the throat, and a linear mean flow distribution Is 
used from the throat to the outlet. If Inlet and outlet critical 
velocity ratios are approximately equal, a linear mean flow 
distribution Is assumed from the Inlet to the throat and from the 
throat to the exit. 

IMN = 3, Does the same thing as IMN = 2, except the outlet flow Is assumed 
to be supersonic, and a sonic flow Is assumed at the throat. 

NMN - Number of values used to describe the flow path. Minimum number 2, 
maximum number 50. 

ISL - Compressibility correction. 

ISL = 0, Physical velocity Is calculated assuming no density variations In 
the tangential direction. 

ISL = 1, A correction Is made In calculating the physical velocity from 

the transformed velocity. This correction Is used to account for 
large variations In density In the tangential direction and gives 
better results for some problems where supercritical flow Is 
encountered . 

Item 3 - Absolute flow conditions at the cascade Inlet station, which Is taken 
to be the leading edge of the farthest upstream body. Any consistent set of 
units may be used. (Format 8F10.0) 

GAMMA AR PTAI TTAI 

GAMMA - Ratio of specific heats. 

AR - Gas constant. 
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PTAI - Total pressure at the Inlet station, measured In the absolute 
reference frame. 

TTAI - Total temperature at the Inlet station, measured In the absolute 
reference frame. 

Item 4 - Inlet relative velocity array. Units should be consistent with 
contents of Item 3. NCS values are to be given. Values are to be taken at 
the blade-row-lnlet (leading-edge) station. (Format 8F10.0) 

WIN(l) . . . WIN(NCS) 

WIN - Inlet relative velocity array. 

Item 5 - Inlet relative flow angle array. Given In degrees. NCS values are 
to be given at the blade-row-lnlet (leading-edge) station. (Format 8F10.0) 

ALPIN(l) . . . ALPIN(NCS) 

ALPIN - Inlet flow angle array. 

Item 6 - Outlet relative flow angle array. Given In degrees. NCS values are 
to be given at the blade-row-outlet (tral ling-edge) station If ICSCO =1. If 
ICSCO = 0, 2, or 3, a dummy line or card must be Included In the 
data set. (Format 8F10.0) 

ALPOUT(l) . . . ALPOUT(NCS) 

ALPOUT - Outlet flow angle array. 

Item 7 - Circulation control array. NBODY values to be given. If ICSCO = 1, 
the fraction of total blade-row circulation to be applied to each element Is 
given. If ICSCD = 3, the tralllng-edge flow angle to be used In the Kutta 

condition Is given In degrees. If ICSCO = 0 or 2, a dummy or blank line Is 

given. (Format 8F10.0) 

CRCLTN(l) . . . CRCLTN(NBODY) 

CRCLTN - Circulation control values for each blade element. 

Item 8 - Blade-row rotation-rate array. Given In revolutions per minute 

(rpm). NCS values to be given. Must be given even If the rotation rate Is 

zero. (Format 8F10.0) 

OMEGA(l) . . . OMEGA(NCS) 

OMEGA - Angular velocity array. 

Item 9 - Flow path description. Radius and stream sheet thickness arrays as a 
function of the meridional coordinate for IMN = 0, 2, and 3. Radius stream 
sheet thickness and mean-flow critical-velocity-ratio arrays as functions of 
the meridional coordinate for IMN = 1. Minimum number 2, maximum number 50. 
(Format 8F10.0) 
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9a. MRSP(l) . . . MRSP(NMN) 

MRSP - Array of meridional coordinates where flow path description values 
are to be given. NMN values are to be given. The far-stream 
boundary-condition stations are taken at MRSP(l) and HRSP(NMN). 

9b. RMSP(l) . . . RHSP(NMN) 

RMSP - Array of flow path transverse radii measured between the 

turbomachine centerline and the blade-to-blade surface of revolution 
at the corresponding MRSP array locations. NMN values to be given. 


9c. BESP(l) . . . BESP(NMN) 


BESP - Array of stream sheet thicknesses corresponding to MRSP array 

locations. Values may be normalized by any arbitrary value. NMN 
values to be given. 


9d . WMSP( 1,1) . . . WMSP(NMN,1 ) 


WMSP(1 ,NCS) 


WMSP(NMN.NCS) 


WMSP - Arrays of mean critical velocity ratio corresponding to the MRSP 
locations and the the flow case being analyzed. Given when 
IMN = 1. NMN * NCS values to be given. 

Body/Blade Shape Input - Items 10 to 12 are read for each body or blade 
element In the flow. The sequence Is repeated NB00Y times. 


Item 10 - Body-shape Input controls. (Format 1615) 

NELM IBLD ITE NBRK 

NELM - Number of panel elements (number of body points minus one). If IBLD 
and ITE = 1, two panels are added to the value of NELM. Maximum 
number 98 or when IBLD and ITE = 1, maximum number 96. 

IBLD - Type-of-blade description. 

IBLD = 0, Descrete points are used to describe the body. 

IBLD = 1, Blade shape generator used to develop blade coordinates. 

\ • \ . 

ITE - Type of trailing edge. 

ITE = 0, Sharp trailing edge. 

ITE = 1, Rounded trailing edge. 

NBRK - Number of surface slope discontinuities. Maximum number 2. A 
surface slope discontinuity occurs where the body shape makes an 
abrupt change In direction (e.g. , the sharp trailing edge of the 
airfoil In fig. 1 ) . 
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Item 11 - Index of body points where surface slope discontinuities occur. The 
Index of a body point Is the number It Is given when counting the points, 
starting at the body trailing edge and moving around the body In a clockwise 
direction. In figure 1, IBRK(l) = 1, Indicating that a surface slope discon- 
tinuity occurs at the first element end point ( 1 . e . , the trailing edge), and 
IBRK(2) = K, Indicating that a surface slope discontinuity also occurs at the 
kth end point ( 1 . e . , the leading edge In the figure). If NBRK = 0, a dummy or 
blank line must be Included. If IBLO = 1 (see Item 12, Option 2 below) and 
the blade has a sharp leading edge, the Index value for the leading edge Is 
calculated and need not be specified. For blade shapes, always run with NBRK 
='l and IBRK(l) = 1, even for rounded trailing edges. This seems to give the 
best results. (Format 1615) 

IBRK(l) . . . IBRK(NBRK) 

IBRK - Array of body point Indices where slope discontinuities occur. 

Item 12 - Blade geometry Input. Option 1 ( IBLD = 0). Discrete blade surface 
coordinates are given by the user. The coordinates must be given starting at 
the lower-surface trailing edge and moving around the surface In a clockwise 
direction. See figure 1. NELM plus one points are to be given with the first 
and last points being coincident. (Format 8F10.0) 

12a. M(1 ) . . . M(NELM + 1) 

M - Array of blade-surface meridional coordinates. 

12b. PH 1(1 ) . . . PHI(NELM + 1) 

PHI - Array of blade-surface tangential coordinates. 

Item 12 - Blade geometry Input. Option 2 (IBLD =1). Blade shape generator 
will be used to develop blade surface coordinates. The blade generator uses 
spline curve fits of the user-supplied blade-surface coordinates to develop 
the upper and lower blade surfaces. These surfaces are then joined to form a 
sharp leading or trailing edge, or are patched together using either a circle 
or an ellipse to make a blunt leading or trailing edge. The blade surface Is 
then divided Into surface elements by the generator according to the user's 
specifications. See figures 2 and 3. (Format 8F10.0) 

12a. MLE PHILE RA RB TAU 

MLE.PHILE - Coordinates of the leading edge for a sharp leading edge. 

Coordinates of the center of the leading edge circle or 
ellipse for a rounded leading edge. 

RA, RB - Semi-major and semi-minor axes for an elliptical leading edge. 

RA = RB for a circular leading edge. RA = RB = 0 for a sharp 
leading edge. RA and RB are determined from a drawing of the 
leading edge on an M, R * PHI plot. See figure 2. 

TAU - Angle, In degrees, of the orientation of the major axis of the 

ellipse with respect to the m axis on an M, R * PHI plot. See 
figure 2. TAU = 0 for circular and sharp leading edges. 
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12b. MTE PHITE RTE 


MTE.PHITE - Coordinates of the trailing edge for a sharp trailing edge. 

Coordinates of the center of the trailing edge circle for a 
rounded trailing edge. 

RTE - Trailing edge circle radius as found on an M, R * PHI plot. 

See figure 2. RTE = 0 for a sharp trailing edge. For blunt 
trailing edges, two elements are located on the trailing edge 
circle. These elements are then added to the number of ele- 
ments given In Item 10 to bring the total element count to 
NELM + 2. For sharp trailing edges, the size of the first 
tralllng-edge element Is set equal to the size of the leading- 
edge element. 

12c. PLE PUPS FLE FTE 

PLE - Fraction of total blade elements (NELM In Item 10) to be used In the 
leading-edge description for a rounded leading edge. Size of the 
leading-edge element, as fraction of chord, for a sharp leading 
edge. See figure 3. 

PUPS - Fraction of total blade elements minus the leading-edge elements to 
be used In the upper surface description. See figure 3. 

FLE - Element multiplication factor for transition of element sizes from 
leading edge to upper and lower surface elements. Value should be 
taken between 1 and 2. See figure 3. 

FTE - Element multiplication factor for transition of element sizes from 
trailing edge to upper and lower surface elements. Value should be 
taken between 1 and 2. See figure 3. 

1 2d . BTUPS BTUPE SPLUP 

BTUPS - Upper-surface tangency angle with leading edge. Given In degrees. 

Angle to be taken from an M, R * PHI plot. See figure 2. 

BTUPE - Upper-surface tangency angle with trailing edge. Given In degrees. 
Angle Is to be taken from an M, R * PHI plot. See figure 2. 

SPLUP - Floating point number of spline points to be given on the upper 

surface. Minimum number 2, maximum number 50. 

1 2e. M( 1 ) . . . M(SPLUP) 

M - Array of upper-surface meridional coordinates to be used In spline 
curve fit. 

1 2f . PHI ( 1 ) . . . PHI(SPLUP) 

PHI - Array of upper-surface tangential coordinates to be used In spline 
curve fit. 
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Note 1 : The first and last spline points are calculated by the code 

so that these values need not be supplied by the user. However, 
blank spaces for these values are required In the spline-point 
coordinate arrays. 

12g. BTLWS BTLWE SPLLW 

BTLWS - Lower-surface tangency angle with leading edge. Given In degrees. 
Angle to be taken from an M, R * PHI plot. See figure 2. 

BTLWE - Lower-surface tangency angle with trailing edge. Given In degrees. 
Angle to be taken from an M, R * PHI plot. See figure 2. 

SPLLW - Floating point number of spline points to be given on lower surface. 
Minimum number 2, maximum number 50. 

1 2h. M( 1 ) . . . M(SPLLW) 

M - Array of lower-surface meridional coordinates to be used In spline 
curve fit. 

121. PHI ( 1 ) . . . PHI(SPLLW) 

PHI - Array of lower-surface tangential coordinates to be used In spline 
curve fit. 

Note 2 : The first and last spline points are calculated by the code 

so that these values need not be supplied by the user. However, 
blank spaces for these values are required In the spline-point 
coordinate arrays. 

Items 10 to 12 are repeated for each additional body or blade shape until 
NBODY shapes have been described. 

Item 13 - Cascade configuration parameters. See figure 4. (Format 8F10.0) 
REFLGT CHORD STGR S 

REFLGT - Reference length Is used to scale the meridional coordinate In 
the solution output listing (l.e., the meridional coordinate Is 
divided by REFLGT In the output listing). 

CHORD - Chord length. 

STGR - Cascade stagger Is measured as the difference between the PHI 

location of the farthest upstream leading edge and the farthest 
downstream trailing edge. Given In radians. See figure 4. 

S - Cascade pitch. Given In radians. The tangential spacing between 

a common point on the Input body geometry and the repeat of the 
Input geometry In cascade. See figure 4. 
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Item 14 - Output control. (Format 1615) 

IOUT IPLOT 

IOUT - Control for printed output. 

IOUT = 0, Gives echo check, error messages, and compressible flow 
results. 

IOUT = 1, Blade geometry Is added to the printout. 

IOUT = 2, Auxiliary flow output Is added to the printout. 

IOUT = 3, Transformed solution Is added to the printout. This 
Information Is primarily used to debug a solution. 

IPLOT - Control for unformated graphical file output. 

IPLOT = 0, No graphics file Is generated. 

IPLOT = 1, Graphics file Is generated and sent to unit 60. 


Basic Input Notes 

Items 1 to 14 given In the previous section make up the basic Input file. 
This data set supplies the Information which Is needed to calculate the flow 
velocity on the surfaces of the bodies In the flow. It represents the minimum 
Input required to run the code. The following notes are Intended to give the 
user additional Information on how to use the code Input to set up special 
flow cases, and to Indicate to the user some possible problems that may occur 
during a solution: 

(1) When the Kutta condition option (ICSCD = 2 or 3) Is used to determine 
the blade row circulation and downstream flow angle, an Iteration procedure Is 
employed. This Iteration procedure has difficulty converging when the down- 
stream relative conditions are near sonic, and will fall to converge when the 
downstream flow conditions are supersonic. If the Iteration count exceeds 
101, the code stops Iterating and continues on with the calculation, using the 
values of the solution variables at that point. The user should, therefore, 
always check the Iteration count to make sure he has a valid solution. The 
Iteration count Is printed In the compressible-flow output when the Kutta con- 
dition option Is used. 

(2) Plane cascade flows may be calculated by setting all the elements In 
the mean flow path radius array (RMSP) to 1.0, and then considering the M and 
PHI coordinates used to describe the blade shapes to be Cartesian coordinates. 

(3) Incompressible flow calculations are set up using Input flow conditions 
(Item 3). The parameters PTAI and TTAI are still set to the absolute Inlet 
total pressure and temperature, but the ratio of specific heats Is used as a 
calculation control switch and Is set to a value greater than 1000. The remain- 
ing parameter In Input Item 3 (the gas constant AR) Is set so that the fluid 
density Is equal to the expression PTAI/(AR * TTAI). All other Input Items 
remain unchanged. For .Incompressible calculations, IOUT In Input Item 14 
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should be set equal to 2. This will cause the dimensional velocity and pressure 
to be printed out. These values are more useful than the standard compressible 
flow output for Incompressible solutions. 

(4) The user may simulate viscous loss and blockage effects by reducing 
thd stream sheet thickness array (BESP) In the mean flow path description. 

The BESP array must be smoothly reduced from the blade-row Inlet to exit by a 
percentage equal to the amount of total pressure loss expected across the 
blade row. The flow velocities calculated will reflect the effect of the 
viscous blockage; however, the pressures listed by the solution will be In 
error. 


Field-Point Input File 

If the user wants to calculate flow velocities In the flow field ( 1 . e - , 
off the body surface), flow-field Input data Is added to the basic Input file 
described previously. The coordinates of the points where the velocity Is to 
be calculated are given along with reference flow angles at the points. The 
off-body Input Is grouped Into blocks of data. The user may add as many of 
these off-body data blocks as needed. The format of the data blocks Is given 
as follows: 

Item lob - Number of off-body (flow-field) point coordinates to be given In 
data block. Maximum number 98. (Format 1615) 

Item 2ob to Nob - Off-body point data. Nob equals the number given In 
Item lob plus one. (Format 8F10.0) 

MOB PHIOB ANG 

MOB - Meridional coordinate of the off-body point where the flow-field 
velocity Is to be calculated. 

PHIOB - Tangential coordinate of the off-body point where the flow-field 
velocity Is to be calculated. (Given In radians.) 

ANG - Reference angle. Given In degrees. The velocity Is calculated along 
and normal to a line through the off-body point. The reference 
angle orients this line with respect to the meridional direction. 

If the reference angle (ANG) was aligned with the flow, then the 
calculated normal velocity would be zero. If ANG = 0, the flow 
angle will be calculated with respect to the meridional direction. 

Items lob to Nob comprise one block of off-body Input data. The user may 
add as many of these blocks as he wants by simply repeating Items lob through 
Nob for each additional block. The number of off-body points can vary from 
block to block; however, the number of points given In any one block can not 
exceed 98. No separate echo check Is made of the off-body Input data, but 
the data Is printed out with the output listing of the calculated off-body 
velocities . 
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Off-Body Input Notes 

Repeated blocks of Items lob to Nob make up the off-body calculation Input 
file. Information given In this data set locates the points In the flow field 
where the user wishes to calculate the flow velocity. The following notes are 
Intended to help the user In making these off-body velocity calculations: 

(1) For multiple-flow-case runs (NCS >1), only the last block of off-body 
points are stored for solutions 2 through NCS. 

(2) A singularity In the off-body velocity calculations occurs when an 
off-body point Is located at a distance of n*p1tch above or below a control 
point located on an Input body (n = +l,+2, .. ., -1,-2, ...). As the off-body 
points approach this location, errors will develop In the velocity calculation. 
These errors can be avoided by restricting the off-body points to a region 
which Is within a tangential location of plus or minus 0.9 * pitch of the 
Input body location as shown In figure 5. 


Sample Input 

A sample Input data set for a three-body problem (two cylinders In a 
channel) Is shown In Table I. Labels for the various sections of the Input 
have been Inserted between the Input data lines. The Input lines that make up 
the actual data set have been appended with Item numbers which correspond with 
the descriptions given previously. 


OUTPUT FILE DESCRIPTION 

The output that Is generated by the multiple-body-panel code Is controlled 
by Item 14 In the basic Input file. Two separate files may be generated. 

These are the listing or printed output file of the tabulated flow solution, 
and an unformatted flow solution file, which Is Intended for use with a post- 
processing graphics program. 


Output Listing File 

The code output listing file Is sent to unit 06, which Is the line printer 
on most computer systems. It Is a formatted data set. If the maximum amount 
of output Is requested by setting I0UT = 3 In the basic Input file (Item 14), 
the following Information appears on the listing In the order given below: 

(1) Labelled echo check of basic Input file. 

(2) Transformed and physical blade shape geometry plus any blade shape 
generator code diagnostics. Listed for each body. 

(3) Code-generated solution diagnostics. 

(4) Transformed velocity solution and singularity strengths for each body. 

(5) Compresslble-f low-solution listing and auxllary-f low-solution output 
for each body. 
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(6) Off-body (flow-field) velocity-calculation 1 results . 

For multiple-flow-case runs (NCS > 1, Item 1 In basic Input file), Items 3 
to 6 are repeated for each additional flow case. The solution generates the 
output listing at various stages In the calculation. This allows the user to 
receive echo check and blade geometry even though the flow solution falls. 


Graphics File 

Solution Information, which can be used In a user-generated graphics code. 
Is sent to unit 60 by the code. The graphics file Is generated If IPLOT = 1 

In Item 14 of the basic Input file. It Is an unformatted data set. The file 

contains the following Information In the order given below: 

(1) The number of bodies In the calculation (NBODY In Item 1 basic Input 

file), the number of flow cases In the run (NCS In Item 1 basic Input file), 

the number of the body with the largest meridional coordinate (bodies are num- 

bered according to the order of their Input In the Input file), the number of 
the body with the smallest meridional coordinate, the mean-f low-calculation 
control parameter (IMN In Item 1 basic Input file), the number of points used 
to describe the mean flow (NMN In Item 1 basic Input file), cascade pitch (S In 
Item 13 basic Input file), and the calculation reference length (REFLGT In 
Item 13 basic Input file). 

(2) The number of elements used to describe each body (NELM In Item 10 
basic Input file) . 

(3) Index of body coordinate with the minimum meridional value for each 
body. Index refers to the the location of the point In the coordinate array 
as sent to this file. 

(4) Index of body coordinate with the maximum meridional value for each 
body. Index refers to the the location of the point In the coordinate array 
as sent to this file. 

(5) The meridional coordinate array of the body control points where the 
surface velocities are calculated. 

(6) The meridional coordinate array of the panel/element end points that 
make up the body. 

(7) The tangential coordinate array In radians of the body control points 
where the surface velocities are calculated. 

(8) The tangential coordinate array In radians of the panel/element end 
points that make up the body. 

(9) Local transverse radius array at the panel end points. 

(10) The panel end point transformed meridional coordinate array. 

Items (5) through (10) are repeated for each body as they are ordered In 
the Input. 
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(11) Meridional coordinate array for mean flow data (MRSP In Item 9 basic 
Input file). 

(12) Transverse radius array at mean-flow meridional coordinates (RMSP In 
Item 9 basic Input file). 

(13) Stream sheet thickness array at mean-flow meridional coordinates 
(BESP In Item 9 basic Input file). 

(14) Mean-flow critical velocity ratio array at the mean flow meridional 
coordinates (Either calculated or given In Item 9 basic Input file). 

(15) Index of the body control-point coordinate which marks the stagnation 
point on the body. Index refers to the the location of the point In the body 
coordinate array as sent to this file. If this Index Is ILU, then control 
points ILU to NELM lie on one side of the stagnation point while control 
points ILU - 1 to 1 lie on the other side. 

(16) Distance along the body surface of the control points from the 
leading-edge stagnation point. 

(17) Local relative critical velocity ratio at the body control points. 

(18) Ratio of local static pressure to local relative total pressure at 
the body control points. 

(19) Pressure coefficient at the body control points. Reference dynamic 
head Is based on the upstream-station relative velocity and density. 

(20) Local relative Mach number at the body control points. 

(21) Relative static pressure at body control points. 

(22) Relative velocity at body control points. 

Items (15) through (22) are repeated for each body. 

(23) Upstream and downstream relative velocity of the mean flow. 

(24) Number of off-body points where the flow field velocity Is to be 
calculated (Item lob In field point Input file). 

(25) The transformed meridional coordinate of the off-body point, the 
tangential coordinate for the off-body point, the relative velocity at the 
off-body point, and the flow angle In radians with respect to the Input 
reference angle. 

For the first flow case, Items (24) and (25) are repeated for all the 
blocks of off-body points given In the field-point Input file. For all 
remaining flow cases, only the last block of off-body points are written out. 
(See field-point Input file description for more Information.) 

(26) The Integer value "9999" Is written to mark the end of one flow case 
and the start of data for the next. 
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Items (14) to (26) are repeated for each flow case. 

The graphics file Is generated by write statements In three subroutines. 
Items (1) to (13) are written out In the subroutine GEQOPT. Items (14) to 
(23) are written out In subroutine CMPOPT. Items (24) to (26) are written out 
In subroutine 0FFB00. If the user plans to use this file, he should look at 
these write statements In order to get a better Idea of how to go about reading 
In the Information In the file. 


TEST CASES 

Three sample test-case Input data sets accompany the code. The first case 
Is an Isolated, two-element airfoil. The calculated surface-velocity distri- 
bution and geometry are given In figure 6, along with an exact solution by 
Williams (ref. 2). The second case Is a centrifugal compressor Impeller with 
splitter blade. Figure 4 shows the problem geometry, and figure 7 shows the 
calculated surface velocity. A flow solution from the finite difference code 
TANSON (ref. 3) Is also given on the figure. The final test case Is a three- 
body problem of two cylinders In a channel. The off-body-velocity calculation 
capability of the code was used to generate the vector plot of the solution 
flow field shown In figure 8. 
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TABLE I. - SAMPLE INPUT FILE. {ONLY LINES APPENDED WITH AN ITEM NUMBER ARE INCLUDED IN THE INPUT.) 


BASIC INPUT FILE 


INPUT DATA SET TITLE 


INCOMPRESSIBLE FLOW THROUGH CYLINDERS - DENSITY 


ITEM 1 


CALCULATION CONTROL PARAMETERS 

311020 ITEM 2 

FLOW CONDITIONS 


1000 . 

2. 2000. 

500. 


ITEM 

3 

10. 




ITEM 

4 

0. 




ITEM 

5 

0 . 




ITEM 

6 

0. 

0.5 0.5 



ITEM 

7 

0. 




ITEM 

8 

MEAN FLOW 

PATH DESCRIPTION 





-5. 

7. 



ITEM 

9A 

1 . 

1 . 




9B 

* 1 . 

1 . 




9C 

FIRST BODY 

GEOMETRY DESCRIPTION ( 

OPTION 2) - 

CHANNEL WALL (FLAT PLATE) 



60 

1 1 1 



ITEM 

10 

1 




ITEM 

11 

-3. 

-1. 0.01 

0.01 

0 . 

ITEM 

12A 

5. 

-1. 0.01 




12B 

.05 

.5 1.75 

1.5 



12C 

0. 

0. 2. 




12D 

0. 

0. 




12E 

0 . 

0. 




12F 

0. 

0. 2. 




12G 

0 . 

0. 




12H 

0. 

0. 




121 

SECOND BODY GEOMETRY DESCRIPTION 

(OPTION 1) 

- CYLINDER 




40 0 0 0 ITEM 10 

(BLANK/DUMMY LINE) ITEM 11 

0.500000 0.493840 0.475530 0.445500 0.404510 0.353550 0.293890 0.227000 ITEM 12A 

0.154510 0.078220 0.000000 -0.078210 -0.154500 -0.226990 -0.293890 -0.353550 12A 

-0.404500 -0.445500 -0.475520 -0.493840 -0.500000 -0.493850 -0.475530 -0.445510 12A 

-0.404520 -0.353570 -0.293910 -0.227010 -0.154530 -0.078240 -0.000020 0.078190 12A 

0.154490 0.226970 0.293870 0.353530 0.404490 0.445490 0.475520 0.493840 12A 

0.500000 12A 

0.000000 -0.078220-0.154510 -0.227000 -0.293890 -0.353550 -0.404510 -0.445500 ITEM 12B 
-0.475530 -0.493840 -0.500000 -0.493840 -0.475530 -0.445510 -0.404510 -0.353560 12B 



-0.293900 -0.227000 -0.154520 -0.078230 -0.000010 

0.293880 0.353540 0.404500 0.445490 0.475520 

0.475540 0.445510 0.404520 0.353570 0.293910 

0.000000 

THIRD BODY GEOMETRY DESCRIPTION (OPTION 1) - CYLINDER 
40 0 0 0 


(BLANK/DUMMY LINE) 


2.500000 

2.493843 

2.475528 

2.445502 

2.404509 

2.154510 

2.078219 

2.000004 

1.921786 

1.845496 

1.595496 

1.554501 

1.524474 

1.506158 

1.500000 

1.595482 

1.646434 

1.706093 

1.772987 

1.845471 

2.154486 

2.500000 

2.226973 

2.293872 

2.353535 

2.404491 

1.000000 

0.921783 

0.845492 

0.773005 

0.706107 

0.524472 

0.506156 

0.500000 

0.506155 

0.524470 

0.706100 

0.772996 

0.845481 

0.921771 

0 .999987 

1.293879 

1.353540 

1.404497 

1.445495 

1.475521 

1 .475535 
1.000000 

1.445514 

1.404524 

1.353572 

1.293915 


CASCADE CONFIGURATION 

1. 1. 0. 3. 

OUTPUT CONTROL 
rv> 
co 

2 1 

FIELD POINT INPUT FILE - BLOCK ONE 


5 

-2. -.5 0. 

- 2 . 0 . 0 . 

-2. .5 0. 

- 2 . 1 . 0 . 

-2. 1.5 0. 

BLOCK TWO 

3 

-1. -.75 0. 

-1. -.5 0. 

-1. -.25 0. 


END OF INPUT DATA FILE 


0.078200 0.154490 
0.493840 0.500000 
0.227020 0.154540 


2.353553 2.293893 
1.773009 1.706113 
1.506153 1.524467 
1.921762 1.999977 
2.445490 2.475518 

0.646447 0.595492 
0.554494 0.595487 
1.078203 1.154494 
1.493840 1.500000 
1.227020 1.154535 


0.226980 

12B 

0.493850 

12B 

0.078250 

12B 


12B 



ITEM 

10 


ITEM 

11 

2.226995 

ITEM 

12A 

1.646452 


12A 

1.554489 


12A 

2.078195 


12A 

2.493838 


12A 

12A 

0.554497 

ITEM 

12B 

0.646441 


12B 

1.226980 


12B 

1.493848 


12B 

1.078246 


12B 

12B 


ITEM 13 


ITEM 14 


ITEM 1 OB 
ITEM 2 OB 
ITEM 3 OB 
ITEM 4 OB 
ITEM 5 OB 
ITEM 6 OB 


ITEM 1 OB 
ITEM 2 OB 
ITEM 3 OB 
ITEM 4 OB 



Tangential coordinate, PHI 



Figure 1. - Blade shape input, option 1, IBLD - 0. 



Figure 2. - Blade shape input, option 2, IBLD * 1. 



O User's input surface points 
■ Leading edge elements, total number equals PLE * NELM 
a Upper surface elements, total number equals PUPS » (1 - PLE) * NELM 
▼ Lower surface elements, total number equals [1 - PUPS * (1 - PLE) - PLE] * NELM 
• Trailing edge elements, total number equals 2 for blunt trailing edge 



controlled by FTE 

Figure 3. - Output from blade generator, option 2, IBLD • 1. 



Figure 4 - Blade-row geometry Input parameters. 




Bodies at 
-1 pitch 


Figure 5. - Off-body velocity point placement. 






Critical velocity ratio, W/W 
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